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Abstract. We compare the theoretical anisotropic upper critical field Hc{Q,T) of 
a quasi-two-dimensional d-wave superconductor with recent H C 2 data for the layered 
organic superconductor n — (BEDT — TTF)2Cu(NCS)2- We find agreement both 
with regard to the angular and the temperature dependence of He- This supports 
the suggestion that the Fulde-Ferrell-Larkin-Ovchinnikov state (FFLO state) exists in 
this material for exactly plane-parallel orientation of the magnetic field. Indications of 
precursor states, occurring for small deviations from the plane-parallel field direction, 
are also pointed out and further measurements for confirming the existence of the 
FFLO state are proposed. 



PACS numbers: 74.25.Ha,74.70.Kn,74.80.Dm 

If the superconducting state in high magnetic fields is limited by paramagnetic 
pair breaking alone, the transition from the homogeneous superconducting state to the 
normal conducting state may proceed either directly, at the Pauli limiting field Hp, or 
via an interposed inhomogeneous superconducting state predicted in 1964 by Fulde and 
Ferrell [1] and by Larkin and Ovchinnikov [2] (FFLO state). The latter state, which 
stabilizes as a consequence of spin polarization, has attracted considerable interest over 
the years, but has, despite of this, not yet been definitely verified experimentally. 

Favourable conditions for observing the FFLO state are found in clean 
superconductors with orbital critical fields much larger than Hp. In practice, it 
seems always necessary to reduce the orbital pair breaking effect by using layered 
superconductors with nearly decoupled planes or extremely thin films (quasi-two- 
dimensional superconductors) and applying the magnetic field in a direction parallel 
to the conducting planes. Several classes of superconducting materials with favourable 
conditions for observing the FFLO state do exist. These include the "classical" 
intercalated transition metal dichalcogenides as well as more exotic materials like High- 
T c compounds and organic superconductors. 

In this Letter we refer to a recent measurement of the upper critical field H& in 
the organic superconductor k — {BEDT — TTF) 2 Cu(NCS)2 [3]. This layered material 
shows strong anisotropy of the superconducting properties with regard to out-of-plane 
directions. In addition, a number of experiments listed in reference [3] are interpreted 
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in terms of a strong in-plane anisotropy, of d-wave type, of the gap parameter. The 
coherence length £j_ perpendicular to the layers of this clean material is smaller than 
the interlayer spacing d and one expects an extreme reduction of orbital pair breaking 
for plane-parallel applied field. The angular dependence of H c2 was measured in Ref [3] 
both with respect to the angle 9 between applied magnetic field and the direction normal 
to the conducting planes and with respect to the azimuthal angle 0, which denotes the 
direction of magnetic field lying within the plane. The results showed, as expected, a 
strong variation of H c2 with 6. On the other hand, no dependence on was observed. 
The maximal value of H c2 at the plane-parallel position = n/2 was of the order of, 
but 50% higher than, the Pauli paramagnetic limit Hp. These facts led the authors of 
reference [3] to propose that their in-plane critical field is the phase boundary between 
the normal-conducting state and the FFLO state of a d-wave superconductor. 

We examined this question by calculating the angular and temperature dependence 
of the theoretical phase boundary H C (Q,T) between the normal-conducting state 
and the superconducting states of a quasi-two-dimensional d-wave superconductor. 
Comparison with the data of reference [3] showed good agreement, supporting the 
hypothesis that in this experiment the phase boundary of a FFLO state (for d-wave 
superconductors) has been observed for the first time. 

We assume that the coupling between the conducting planes of k — (BEDT — 
TTF) 2 Cu(NCS)2 can be neglected, making our problem effectively two-dimensional. 
Then, if the field has both a perpendicular and a parallel component, the 
superconducting state is limited by both orbital and paramagnetic pair breaking. For 
exactly plane-parallel magnetic field the FFLO phase should be realized below a reduced 
temperature (t = T/T c ) of t w 0.4 [4]. Such a situation, with competition between both 
pair breaking effects, has been studied first by Bulaevskii [5]. His treatment was later 
generalized to arbitrary temperatures and d-wave superconductors by Shimahara and 
Rainer [6]. Below the stability limit of the normal-conducting state, which will be 
referred to as Hc{Q,T), a series of different inhomogeneous superconducting states, 
depending on O, appear. For s-wave superconductors each one of these states belongs 
to a particular value of Landau's quantum number n, which takes integer values 
n — 0, 1, 2, ... These states are (for s-wave) the following: 

• The vortex state for small O belongs to n — 0. 

• A series of inhomogeneous states for O near 7r/2, each one characterized by a single 
value n > 0, with n increasing with increasing O. 

• The FFLO state for O = ir/2, which may be characterized by n — > oo. 

The structure of the higher Landau level states, for n > 0, has been calculated recently 
for s-wave superconductors by minimizing the quasiclassical free energy [7]. For d-wave 
superconductors [6] a state below Hc(Q,T) is no longer characterized by a single value 
of n but rather by an infinite subset {no, no ± 4, no ± 8, . . .}. However, the dominant 
contribution may still be characterized by a single number n, which increases again 
with increasing O and approaches infinity in the FFLO limit. Thus, basically the 
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above classification scheme remains valid for d-wave symmetry. The phase boundary 
of the 'pure' FFLO state for d-wave superconductors [the curve Hc{tt/2,T)] has been 
calculated by Maki and Won [4]. 

The linearized gap equation to be solved is given by [5,6] 

- log(£)A(f)) = nk B T r , J" - d -f [m 2 ]] [1 

T c Jo smh{7!- kbTs) Jo 2tt 1 J 

-cos[ S {^-^n}]]A(f). (1) 

We consider a cylindrical Fermi surface, appropriate for the present two-dimensional 
problem, with the Fermi velocity vp = vp(e x cos 4>+e y sin 0). The gap parameter is given 
by A(r,p) = A(r) r y(p), where 7(p) = 1 for s-wave, and 7(j3) = V2(p x —p 2 y ) = l + cos(40) 
for d-wave pairing. The canonical momentum is defined by II = -(V — ij?-A). 
The magnetic field H is assumed to lie in the yz-plane, with H y = Hn = H sin 
and H z = H± = HcosQ. We use the following gauge for the vector potential: 
A x = Hnz — H±y, and A y = A z = 0. Paramagnetic pair breaking enters via the term 
HoH in (1); the electron's magnetic moment is /j,q = —gL^B/2, with the Lande factor 
~ 2 and Bohr's magneton jiB = he/(2mc). The method of reducing equation (1) to 
a set of algebraic equations follows exactly reference [6] and need not be repeated here. 
As a first check of our numerical method we compared our results with reference [6] and 
found complete agreement whenever the same set of input parameters was used. 

Let us proceed to a comparison of the solutions He of equation (1) with the H c2 
data reported in reference [3] . The first point to address is the independence of H c2 on 
the azimuthal angle <fi. Such a dependence may easily be incorporated in the present 
model by replacing tfp in (1) by vf [e x cos(0' + 0) + e y sin(0' + 0)]. As long as orbital 
pair-breaking is present for a parallel field component, the symmetry-breaking term 
A x = H\\z has to be kept in (1). If, on the other hand, complete decoupling of the 
(infinitely thin) conducting planes can be assumed, the term H\\z can be dropped and 
He becomes independent of 0, as can be explicitely confirmed numerically. In this 
context, we recall that a three-dimensional d-wave superconductor shows anisotropy of 
the upper critical field [8] . Thus, the observed independence of H c2 on is in agreement 
with the present model and, in particular, with the assumption of quasi-two-dimensional 
superconductivity. 

If T is measured in units of T c and H in units of /iA (where A is the BCS- 
gap at T = 0), the present model requires only one single parameter ksTc/Ep to be 
fitted. This parameter is proportional to the "bulk" ratio of the orbital and spin critical 
fields %c/(2e^l) and k B T c / fi respectively. In the present anisotropic model the actual 
ratio of spin and orbital pair-breaking depends on 6 and may be written in the form 
kBT c /(Ep cos ©). The best fit to the G-dependence of H c2 at T = 1.45 K (see figure 
4 of reference [3]) has been obtained for ksTc/Ep = 0.058. This value of ksTc/Ep is 
consistent with a critical temperature T c = 10.4 K of the sample studied in reference 
[3], and a Fermi energy of the order of 100 K as estimated from several experiments [9]. 
Using this value we found very good agreement between H c and the data of reference [3], 
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as shown in figure 1. The theoretical curves in figure 1 have been calculated assuming 
d-wave symmetry; the difference between d-wave and s-wave was found to be small 
except very close to T = K and = 90 deg. 
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Figure 1: Square of upper critical field, normalized to its value at = 0, as a function of 
0. Full squares: data of reference [3] at t — T/T c = 0.14. Full line: theoretical result for 
the Bulaevskii-Shimahara-Rainer phase boundary for t = 0.14. Dashed line: theoretical 
result for t = 0.8. 



The H C 2 data show a small but clearly visible kink near the plane-parallel 
orientation, at ~ 87 deg. A similar feature is also found in the theoretical phase 
boundary Hc{&), as shown in figure 1 (where the square of the critical fields has been 
plotted in order to make this discontinuous change in slope better visible). This kink 
indicates the transition from the vortex state, with n — 0, to the first of the above 
mentioned FFLO-precursor states, with n = 1. Still closer to 90 deg equation (1) 
yields additional transitions corresponding to n — 2, 3, which are still visible in figure 1 in 
the theoretical curve but not in the data points. The H c -curves describing the n — 0, 1, 2 
transitions, for the same material but higher T, are shown on a larger scale in figure 
2. The order parameter structure of the precursor phases, where pairing takes place in 
Landau levels n > 0, has been investigated recently for s-wave superconductors [7]. Two 
types of such precursor states have been found: (i) quasi-one-dimensional states, which 
may be considered as a mixture of rows of vortices and one-dimensional FFLO-type 
oscillations, and (ii) two-dimensional lattices with several zeros of the order parameter 
with different vorticity [7]. Such unusual states [of type (ii)] have been predicted to 
occur in the extremely high field region where quantization of single electron levels 
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becomes important [10]. The present arrangement might provide a relatively feasible 
way of observing such vortex structures. It should be mentioned, however, that for 
d-wave superconductors neither the equilibrium structure of the FFLO state nor that 
of the n > states has been calculated so far. 



k-(BEDT-TTF) 2 Cu(NCS) 2 k B T c /E F =0.058, t=0.4 




76 - 



74 I i i i i i j i 

89.6 89.4 89.2 89 88.8 88.6 

e [deg] 

Figure 2: The branches Uq = 0, 1, 2 of the upper critical field, as calculated from equation 
(1), are plotted in more detail, using the same value of kBT c /E F as in figure 1 but a 
higher temperature t = 0.4. 



The shape of the H C (Q) curve depends distinctively on temperature, as shown 
by the plot (dashed line) of Hc{Q) 2 at t = 0.8 in figure 1. The reason for this 
enhancement at higher T is, of course, that paramagnetic pair breaking becomes less 
effective at higher temperature. Data at higher T have not been reported in reference 
[3] but would be useful in order to check the present interpretation. Looking for similar 
measurements we found old data by Morris and Coleman [11] for intercalated transition 
metal dichalcogenide TaS2 — (pyridine) samples. In this material, which represents a 
nearly perfect realization of two-dimensional superconductivity, an unexplained anomaly 
with regard to the behavior of the upper critical field at different T has been reported 
(figure 10 of reference [11]). We find excellent agreement (see figure 3) comparing the 
data of reference [11] with the solutions of equation (1) for an s-wave superconductor. 
Again, a single parameter has been adjusted (ksTc/Ep = 0.024) to obtain both of the 
theoretical curves shown in figure 3. The resistance data reported in reference [11] 
show a non-monotonic behavior near the plane-parallel field orientation (see figure 3 of 
reference [11]), which may be due to transitions to the n > states. The latter states 
are discussed in more detail in reference [7]. 



6 



TaS 2 -(pyridine) k B T c /E F =0.024 




90 88 86 84 82 80 

e [deg] 



Figure 3: Comparison of the O-dependence of the upper critical field of TaSi — 
{pyridine), as reported in reference [11], with the solutions of equation (1) for s-wave 
superconductivity at t — 2.86 and t = 1.4. 
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Figure 4: Comparison of the temperature dependence of the plane-parallel upper critical 
field reported in reference [3] with the theoretical result reported in reference [4]. 
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Finally, let us compare the measured temperature-dependence of H c2 for the plane- 
parallel field orientation (figure 3 of reference [3]) with Hc{tt/2,T). According to our 
interpretation of these data, the states below H C 2(it/2,T) should be a d-wave version 
of the FFLO state for T < T* = 0.4 T c , and the homogeneous superconducting state 
for T > T*. This phase boundary has been calculated first by Maki and Won [4]. 
The Shimahara Rainer d-wave phase boundary must agree with reference [4] in the 
limit 6 — > 7r/2. We found agreement, except for the steep rise of H c (tt/2,T) below 
0.05T C reported in reference [4] (a possible reason for this discrepancy may be slow 
convergence of our numerical method for low T and high n). The comparison between 
theory [4] and experiment [3] depicted in figure 4 shows again fairly good agreement. A 
characteristic difference in temperature variation above and below T* is visible in the 
data points, although it is less pronounced than in the theoretical curve. The difference 
in critical fields between s-wave and d-wave is again rather small; at T = 0.05 T c , the 
lowest temperature, where measurements for k — (BEDT — TTF) 2 Cu(NC 'S) 2 have 
been reported [3], both are approximately given by the standard 2Z)-result [12] for the 
FFLO state, HqH c = A , which exceeds the Pauli limiting field by = 40%. Thus, 
these numbers do also fit well into the FFLO interpretation of the phase boundary for 
plane-parallel field orientation. 

Summarizing, the proposal of Nam et al. [3] that upper critical field data for a 
plane-parallel field orientation in the layered organic superconductor k — (BEDT — 
TT F) 2 Cu(N C S) 2 should be interpreted in terms of a FFLO state, has been supported 
by our calculations. The data agree with the predictions of a model of a quasi-two- 
dimensional superconductor both with regard to the angular and the temperature 
dependence of the critical field. Further confirmation of this interpretation could 
be obtained by means of measurements at higher temperatures, where paramagnetic 
pair-breaking is strongly reduced. If this interpretation is correct, precursor states 
with interesting properties should appear for applied fields close to the plane-parallel 
orientation. 
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